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We argue that the centrosymmetric C2/c symmetry in BiMnOs is spontaneously broken by antiferromag- 
netic (AFM) interactions existing in the system. The true symmetry is expected to be Cc, which is compatible 
with the noncollinear magnetic ground state, where the ferromagnetic order along one crystallographic axis 
coexists with the the hidden AFM order and related to it ferroelectric polarization along two other axes. The 
C2/c symmetry can be restored by the magnetic field i? 35 Tesla, which switches off the ferroelectric 
polarization. Our analysis is based on the solution of the low-energy model constructed for the 3d-bands of 
BiMnOs, where all the parameters have been derived from the first-principles calculations. Test calculations 
for isostructural BiCrOs reveal an excellent agreement with experimental data. 
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BiMnOs is regarded as one of the prominent multi- 
,ferroic materials, where the ferromagnetic (FM) magne- 
tization is coupled to the ferroelectric polarization, thus, 
giving a possibility to control the magnetic properties by 
applying the electric field and vise versa. The ferromag- 
.nctism of BiMnOa is well established today: the Curie 
|temperature is about 100 K and the largest reported 
magnetization is 3.92 fig per formula unit [I], which is 
Iclose to 4 /is expected for the fully saturated FM state. 
'The ferroelectric hysteresis loop was also observed in 
polycrystalline and thin film samples of BiMnOs [2], 
although the measured polarization was small (about 
0.043 /xC/cm^). The multiferroic behavior of BiMnOa is 
|typically attributed to the existence of two different sub- 
lattices: the stereochemical activity of Bi(6s^) lone pairs 
lis believed to be the origin of the structural distortions 
(having the C2 symmetry) and the inversion symmetry 
breaking, which gives rise to the ferroelectricity, while 
the Mn-sublattice is responsible for the magnetism 
.This point of view dominated over several years and was 
documented in many review articles [4j. 

Nevertheless, in 2007 Belik et al. have reexamined 
the crystal structure of BiMnOa and argued that below 
770 K, it is best described by the space group C2/c, 
which has the inversion symmetry [1]. This finding 
was later confirmed by Montanari et al. [S]. According 
to this new crystal structure information, BiMnOa can 
be only in an antiferroelectric state [6] ~ the situation, 
which is rather common for many distorted perovskites 
and hardly interesting from the practical point of view. 
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Of course, these experimental works raised many ques- 
tions. Particularly, what is the origin of the inversion 
symmetry breaking (if any) and ferroelectric response 
in BiMnOa? 



In our previous work, we have constructed an effec- 
tive Hubbard-type model for the 3d-bands of BiMnOa 
located near the Fermi level [7] . All parameters of this 
model were derived in an ab initio fashion on the basis 
of first-principles electronic structure calculations and, 
apart from approximations inherent to the construction 
of the low-energy model [HIE], no adjustable parame- 
ters have been used. Then, from the solution of this 
model in the mean-field Hartree-Fock (HF) approxima- 
tion we derived parameters of interatomic magnetic in- 
teractions and argued that, in the C2/c structure, the 
nearest-neighbor FM interactions (Jnn) compete with 
the longer range antiferromagnetic (AFM) interactions 
{Jlr, see Fig[T]). The latter make the sites 1 and 2 
inequivalcnt and tend to break the inversion center lo- 
cated in the midpoint of the cube diagonal connecting 
these two sites. Both interactions are directly related to 
the orbital ordering realized in BiMnOa below 550 K. 
This suggests that the ferroelectric behavior of BiMnOa 
could be related to some hidden AFM order, which is 
driven by Jlr and which in [7| was denoted as the tiiT 
order, referring to the directions of spins at four Mn- 
atoms in the primitive cell. But, now we face a more 
fundamental problem: how to reconcile the AFM struc- 
ture, which is needed in order to break the inversion 
symmetry, with the FM behavior of BiMnOa, which is 
clearly seen in the experiment? 
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Fig.l. (a) Crystal structure of BiMnOg. Bi-, Mn-, 
and 0-atoms are indicated by big, medium, and small 
spheres, respectively. Primitive translations are shown 
by arrows. Four Mn-atoms, composing the primitive 
cell of BiMnOs, are indicated by numbers, (b) Distri- 
bution of charge densities associated with the occupied 
eg orbitals (the orbital ordering) realized in the low- 
temperature phase of BiMnOa and related to it nearest- 
neighbor ferromagnetic (Jjvjv) and longer range anti- 
ferromagnetic (Jlr) interactions. The position of the 
inversion center in the C2/c structure is marked by the 
cross. 



The clue to this problem may be related to the quasi- 
degeneracy of the FM and tiit AFM states. Indeed, 
according to the HF calculations, the total energy dif- 
ference between these two states was only 0.5 mcV per 
one formula unit [7j . Although the value itself maybe at 
the edge of accuracy of the model calculations, it gives 
a clear ideas about characteristic energy scale we have 
to deal with in the process of analysis of the magnetic 
properties of BiMnOa. 

In the present work we further develop this idea and 
argue that the AFM tiit a-nd FM states can be eas- 
ily mixed by the rclativistic spin-orbit (SO) interaction 
^§•1. Although the expectation value of the orbital mag- 
netic moment /i^l is small in manganites (typically, few 
hundredths of [10] )j the constant of the SO interac- 
tion ^ is about 50 meV. Thus, the energy gain caused 
by the SO interaction can become comparable with the 
total energy difference between the FM and tiit AFM 
states. First, using general symmetry considerations, wc 
argue that the true magnetic ground state of BiMnOs 
is noncoUinear and includes both FM and AFM compo- 
nents of the magnetic moments. This naturally resolves 
the puzzle of multifcrroicity of BiMnOs: the inversion 
symmetry is broken by the hidden AFM order, while 



the FM component appears due to the canting of the 
magnetic moments. The idea will be further supported 
by the direct optimization of the magnetic structure of 
BiMnOa in the model HF calculations with the SO in- 
teraction. Finally, we show how the ferroelectric po- 
larization can be controlled by the external magnetic 
field, which acts on the FM component of the magnetic 
moments. Particularly, we predict the existence of the 
critical field Be, which restores the inversion symmetry 
and switches off the ferroelectric polarization. 

We begin with the symmetry considerations. The 
space group C2/c has four symmetry operations: Si=E, 
S2=i, ^3={my|a3/2}, 54={C2|a3/2}, where E is the 
unity, / is the inversion, rhy is the mirror reflection of the 
y-axis, and Cy is the 180° rotation around the y-axis. 
The last two operations are combined with the transla- 
tion by the vector 3.^/2 (see FigH^). The directions of 
X, y, and z are related to the directions of the primitive 
translations as follows: a;||(ai+a2+0.34a3), j/||(a2— ai), 
and -z||a3. Four Mn-atoms are divided in two classes: 
(1,2) and (3,4). Atoms of one class are transformed by 
the symmetry operations only to each other, but not 
to atoms of another class. The magnetic groups corre- 
sponding to the C2/c space group are obtained by com- 
bining the symmetry operations S2-Si with the time 
inversion T, which additionally flips the directions of 
the magnetic moments. Then, there arc two magnetic 
groups, which can be regarded as candidates for the 
magnetic ground state of BiMnOs: Gi={Si, S2, S3, S4} 
and G2={5'i, 5*2, ^sT, 5'4T}. The third possibility in- 
cludes the symmetry operation S2T=IT . However, 
since / transforms the atoms 3 and 4 to themselves, 
the symmetry operation IT corresponds to the non- 
magnetic sublattice (3,4). Such a situation, although 
formally possible, cannot be realized as the magnetic 
ground state of BiMnOs, because it contradicts to the 
first Hund rule [7]. Then, x and z projections of the 
magnetic moments obey the same transformation rules 
under the symmetry operations Si-S^. Therefore, the 
magnetic structure of BiMnOs can be generally abbre- 
viated as AB-CV, where each capital letter stands for 
the magnetic arrangement (ferromagnetic - F, antifer- 
romagnetic - A, or with zero projection of the magnetic 
moment - Z) formed by the x and y projections of the 
magnetic moments of the atoms (1,2) and (3,4): first 
two capital letters describe the magnetic arrangement in 
the pair (1,2) formed by the x (A) and y (B) projections 
of the magnetic moments; while second two capital let- 
ters describe the magnetic arrangement in the pair (3,4) 
formed by the x (C) and y (2?) projections of the mag- 
netic moments. For example, the magnetic structure 
corresponding to Gi and G2 can be denoted as ZF-AF 
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and FZ-FA, respectively. Moreover, we consider the so- 
lutions with the spontaneous symmetry breaking corre- 
sponding to three possible subgroups of the space group 
C2/c. From this point of view, the most promising are 
two magnetic subgroups: G3={S'i, 6*3} (the space group 
Cc, No. 9 in International Tables) and G4={Si, S3T}, 
which preserve the atomic classes (1,2) and (3,4). They 
correspond to the magnetic structures AF-AF and FA- 
FA, respectively. In the other words, contrary to Gi 
and G2, the magnetic subgroups G3 and G4 allow for 
the nontrivial AFM arrangement in the pair (1,2) with 
finite projections of the magnetic moments. Such a situ- 
ation becomes possible due to the breaking of the inver- 
sion symmetry / (relative to the center, which is shown 
by the cross in Figdjs) . In the case of the coUincar mag- 
netic arrangement, the magnetic structures AF-AF and 
FA-FA are reduced to either FM or AFM tiit ones, 
which have been considered in [7] without the SO inter- 
action. Other magnetic subgroups, generated by 5*2=/ 
and S'4={C'^|a3/2} destroy the atomic classes and make 
the atoms forming the pairs (1,2) and (3,4) inequiva- 
lent. Although such solutions arc formally possible, all 
of them arc unstable and not realized as the magnetic 
ground state of BiMn03, as it will become clear from 
the solution of the low-energy model. 

Our next goal is the search for the true (and so 
far unknown) magnetic ground state of BiMn03 based 
on the solution of the low-energy model. However, 
before doing any predictions for BiMn03, we would 
like to show some test calculations for the isostruc- 
tural compound BiCr03 [11] , which give some idea 
about the accuracy of our approach. The formal con- 
figuration of Cr'^+ ions in BiCr03 is t^g- Since the 
i2g-levels are well separated form the eg-ones by the 
crystal-field (CF) effects, the orbital degrees of freedom 
in BiCr03 arc quenched and the Hubbard model can 
be further mapped onto the Heiscnberg model Tis ~ 
~Ej>j Jij^i^j with the spin 3/2.^' The parameters 
of such a model are shown in Fig[2^. Generally, the 
distribution of the magnetic moments associated with 
the existence of two inequivalent sublattices (1,2) and 

'''We used the same strategy as in our previous work on 
BiMnOs [?]■ First, we construct an effective Hubbard- type model. 
The main difference of the model parameters from the ones re- 
ported for BiMnOs [7] is the following: the CF-spIitting between 
two eg-levels is small (134 and 185 meV for the sites 1 and 3, 
respectively), while the t2g-eg splitting is large (about 2.1 eV). 
Moreover, the on-site Coulomb repulsion U is larger in BiCrOa: 
2.63 and 2.73 cV for the sites 1 and 3, resectively. Then, we cal- 
culate Jij using infinitesimal rotations near the tTii (or G-typc) 
AFM ground state. In BiCrOs, other magnetic configurations 
yield practically the same set of the exchange parameters Jij, 
meaning that the mapping onto the Heisenberg model is universe 
and does not depend on the magnetic state. 



(3,4) should be nonuniform. The magnitude of this ef- 
fect can be estimated in the mean-field approximation. 
However, it appears that from the magnetic point of 
view, the sublattices (1,2) and (3,4) are nearly equiva- 
lent. Indeed, the behavior of local magnetization at the 
sites 1 and 3 (FiglJla) is practically indistinguishable in 
the whole temperature range below the Neel temper- 
ature Tn (about 204 K in the mean-field approxima- 
tion). Then, a proper (and much better) estimate for 
Tjv can be obtained by assuming that the main exci- 
tations in BiCr03 can be described by (uniform) spin 
waves and employing the rcnormalizcd spin-wave the- 
ory, which takes into account the spacial correlations 
between the spins [T2|. It yields 2^=123 K, which is 
very close to the experimental value r/v=109 K [T^ . 




Fig. 2. (a) Distribution of occupied t29-orbitals at the 
Cr-sites and corresponding parameters of Heisenberg 
model {JijS^, where S=3/2) associated with different 
bonds in BiCrOs. Bi- and O-atoms are shown by big 
and small spheres, respectively, (b) Corresponding tem- 
perature dependence of the local magnetization at the 
sites 1 and 3 in the mean-field approximation, which 
demonstrates that the two curves are practically indis- 
tinguishable, despite the fact that the sites 1 and 3 are 
inequivalent and belong to difi'erent atomic classes. 

Being encouraged by these results, we now turn to 
the problem of the magnetic ground state and the ori- 
gin of the ferroelectric behavior of BiMn03. First, we 
perform the delicate optimization of the magnetic struc- 
ture of BiMn03 based on the HF solution of the low- 
energy model with the SO interaction. The parameters 
of the low-energy model were reported in [7]. The so- 
lution of the low-energy model with the SO interaction 
was performed along the same line as for other distorted 
perovskite oxides considered in the review article [5] . A 
typical iterative procedure is shown in Figl3l In this 
case, we start with the FM solution, where all magnetic 
moments were aligned along the y-axis, switch on the 
SO interaction, and monitor the behavior of the FM mo- 
ments {My=My) as well as two AFM order parameters 
L,^={Ml-Ml)/2 and L^={M^-M^)/2, which were de- 
veloped at the sites 1 and 2 in the process of iterative 
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solution of the HF equations. In these notations, Af* 
is the projection of the magnetic moment at the site 
i parallel to the a-axis, where a= x, y, or z. In the 



magnetic moments 'M.= {Mx, My, M^) at the sites 1-4 
(in i^lb): 
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Fig.3. Convergence of the ferromagnetic moment 
My = My = My, two antiferromagnetic order parame- 
ters = (Mi-M^)/2 and = (Mi-Mf)/2, and 
the total energy in the Hartree-Fock calculations af- 
ter switching on the spin-orbit interaction. The ini- 
tial configuration corresponded to the self-consistent 
Hartree-Fock solution for the ferromagnetic state with- 
out spin-orbit interaction, where all magnetic moments 
were aligned along the y-axis. The shaded area is the 
region corresponding to the change of the symmetry 
Gi Gg of the Hartree-Fock solution by the antifer- 
romagnetic interactions. 

process of first approximately 10^ iterations, the solu- 
tion obeys the Gi symmetry, and the iterative proce- 
dure is accompanied by the small but steady decrease 
of the total energy (following the sharp drop right after 
switching on the SO interaction). As expected, the x- 
and z-projections of the magnetic moments at the sites 
1 and 2 are exactly equal to zero in this region (so as 
the AFM order parameters Lx and L^)- However, in 
the process of next 10^ iterations we clearly observe the 
lowering of the magnetic symmetry Gi =J> G3, which is 
accompanied by the growth of AFM order parameters 
Lx and (at the expense of the FM magnetization 
along y) and the step decrease of the total energy of the 
system (by about 10 meV per one formula unit). This 
solution finally converges to the following values of the 
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We also tried to start with other magnetic configu- 
rations (including all magnetic configurations consid- 
ered in |7j with different directions of the magnetic mo- 
ments). However, all of them finally converged to the 
solution described above. For example, by starting with 
the AFM tiit configuration and aligning all magnetic 
moments along y-axis, one can obtain a metastable so- 
lution having the G2 symmetry, which is only 3 meV 
higher that the G3 solution discussed above. However, 
even this metastable solution finally converged to the 
G3 one, which seems to be the true global minimum 
of the total energy of the system. These calculations 
clearly show the advantages of working with the low- 
energy models, because in order to reach the true global 
minimum with the SO interaction, we typically need a 
huge number of iterations, which is only affordable in 
the model approach. 

We would also like to emphasize that the sym- 
metry of the magnetic ground state, obtained in the 
present work, is incompatible with the noncentrosym- 
metric crystal structure C2, which was proposed in ear- 
lier studies [14] and was shown to be unstable [6] with 
respect to the centrosymmetric C2/c structure [T]. In- 
stead, we propose that, at least below the magnetic 
transition temperature, the crystal structure of BiMnOa 
should have the Cc symmetry. This finding should be 
checked experimentally and we hope that our work will 
stimulate further activities in this direction. 

The magnetic group G3 allows for the net electric 
polarization P in the plane zx, which is a normal vec- 
tor and does not changes its sign under S'3={my|a3/2}. 
On the other hand, M is an axial vector and both AFM 
order parameters and Lz change the sign under S3. 
Therefore, the ferroelectric polarization caused by the 
magnetic degrees of freedom must be at least bilinear 
with respect to and Lz |15| : 



Pa — ^ XabcLbLc, 



(1) 



where the symbols a, b, and c denote the z and x pro- 
jections of the vectors. Since L^ and Lz are related to 
the the FM magnetization My (due to the conservation 
of |M|), the value of P can be controlled by the ex- 
ternal magnetic field B=(0, i?,0), applied along y and 
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coupled to My. FigH] shows results of HF calculations 
in the external magnetic field. The corresponding in- 
teraction term is described by Hb = — MsB • (2s+l). 
As expected, B saturates the FM magnetization My. 
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operation (the mirror reflection of the y-axis with sub- 
sequent translation by 3.^/2), which allows for the ex- 
istence of both ferromagnetic order (along the y-axis) 
and ferroelectric polarization (in the plane perpendicu- 
lar to y). The full C2/c symmetry can be restored in 
the external magnetic field, which can be used in prac- 
tical applications in order to completely switch off the 
ferroelectric polarization. 

The work of IVS is partly supported by Grant-in- 
Aid for Scientific Research in Priority Area "Anomalous 
Quantum Materials" and Grant-in- Aid for Scientific Re- 
search (C) No. 20540337 from the Ministry of Educa- 
tion, Culture, Sport, Science and Technology of Japan. 
The work of ZVP is partly supported by Dynasty Foun- 
dation, Grants of President of Russia MK-3227.2008.2, 
and scientific school grant SS- 1929. 2008. 2. 



Fig. 4. Beliavior of the ferromagnetic magnetization 
My = My = My and the antiferromagnetic order pa- 
rameters = {Ml-Ml)/2 and = (Mi-Mf)/2 in 
the external magnetic field along the i/-axis. 

Then, other two projections of the magnetic moments 
will decrease, so as the AFM order parameters and 
Lz. Therefore, the ferroelectric polarization ([1]) should 
also decrease. When the field exceeds some critical value 
Be ^ 35 Tesla, both and vanish. Thus, B > Be 
restores the full symmetry of the system, G3 =J> Gi, and 
switches off the ferroelectric polarization. This find- 
ing is very important for practical applications as it 
clearly demonstrates how the ferroelectric behavior of 
BiMnOs can be controlled by the magnetic field. Of 
course, the obtained Be is too large for any practical 
applications. Nevertheless, this theoretical value may 
be overestimated because the model itself does not in- 
clude several important ingredients such as the mag- 
netic polarization of the oxygen sites [7]. The latter 
will additionally stabilize the FM interactions [71 [T6] 
and therefore decrease Be. On the other hand, the FM 
interactions at the oxygen sites will be partly compen- 
sated by AFM interactions, which arc additionally sta- 
bilized by correlation effects beyond the HF approxima- 
tion [16]. Moreover, the relative strength of FM and 
AFM interactions (and therefore Be) could be tuned by 
the external factors such as pressure, defects or the lat- 
tice mismatch in the thin films of BiMnOa deposed on 
different substrates [7]. 

In conclusion, we predict that the centrosymmetric 
C2/c symmetry in BiMnOs is spontaneously broken by 
hidden AFM interactions, which downgrade the actual 
symmetry to the Cc one (No. 9 in International Tables). 
The space group Cc has only one nontrivial symmetry 
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